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We present a microscopic model of interacting fermions where the ground state degeneracy is 
topologically protected. The model is based on a double-wire setup with local interactions in a 
particle number conserving setting. A compelling property of this model is the exact solvability 
for its ground states and low energy excitations. We demonstrate the appearance of topologically 
protected edge states and derive their braiding properties on a microscopic level. We find the 
non-abelian statistics of Ising anyons, which can be interpreted as Majorana-like edge states. 
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Topologically protected ground state degeneracies in 
many-body quantum systems, and the closely related 
(non-abelian) anyonic statistics, are of special interest 
from a theoretical point of view [1, 2], and have been rec¬ 
ognized as promising concepts for scalable fault-tolerant 
quantum computation [3, 4]. A well understood class 
are topological states with Majorana zero-energy edge 
modes appearing within mean-field descriptions of topo¬ 
logical superconductors [5]. These free fermion theories 
have been classified exhaustively [6, 7], and the properties 
of the Majorana zero modes at boundaries [8] and in vor¬ 
tices [9] have been characterized. In contrast, interacting 
and gapless phases are less well understood [10-13], and 
to which extent existence and non-abelian properties of 
edge states carry over to interacting theories is an inter¬ 
esting question lacking conclusive answers [14]. 

The understanding of topological states is driven by 
exactly solvable microscopic models; the paradigmatic 
one for the existence of topologically protected Majo¬ 
rana edge modes is the one-dimensional Majorana chain 
[8]. It has inspired a variety of proposals for its exper¬ 
imental realization in condensed matter systems [5, 15- 
17], and signatures consistent with Majorana modes 
have been experimentally observed [18-21]. Neverthe¬ 
less, these models require large reservoirs to justify their 
mean-held description, whereas very little is known about 
the fate of Majorana zero-energy edge modes in intrinsi¬ 
cally interacting and particle conserving settings. Pre¬ 
vious attempts for number-conserving theories featur¬ 
ing Majorana-like edge states relied either on bosoniza- 
tion [22-26] or on numerical methods (DMRG) [27], while 
the only exactly solvable models require unphysical long- 
range interactions [28]. 

In this letter, we present a microscopic, number- 
conserving theory with local interactions that features 
non-abelian edge states at the boundaries. The theory 
allows for an exact derivation of its many-body ground 
state as well as its low energy excitations, and thereby 
provides a viable playground for analyzing its character¬ 
istic properties. We hnd that the ground state is charac¬ 
terized by a condensate of p-wave pairs with a topolog¬ 
ical degeneracy. The Green’s function exhibits a revival 



Figure 1. Setup, (a) We consider a double chain (two-leg 
ladder) of spinless fermions with upper/lower chain denoted 
as a/b. (b) The number-conserving Hamiltonian is given by 
intra-chain terms Af (x = a, b) and inter-chain couplings Bj. 

at the edges, indicating the appearance of edge states. 
Remarkably, the system can be extended to arbitrary 
wire networks, which allows us to derive the non-abelian 
braiding statistics of the edge states on a microscopic 
level. 

We consider a double chain (two-leg ladder) of spinless 
fermions with L lattice sites. The fermionic creation op¬ 
erators at site i are described by a| (upper chain) and 
(lower chain), see Fig. 1. The many-body Hamiltonian 
H = -I- + 77“*' describing the interacting fermion 

theory combines intra-chain contributions (x = a, b) 
as well as interactions 77“*' between the two chains. The 
intra-chain Hamiltonian takes the form 

=Eti + ( 1 ) 

with the single-particle hopping terms 

Aj = -l- Oi-i-iaJ, Aj = . (2) 

Gonsequently, it combines single particle hopping with a 
nearest-neighbor attraction nf + — The 

inter-chain interaction 77“** takes a similar form 

= -^^(1 + 5.) (3) 

with the pair-hopping between the two chains 
Bi = alal^.^bibi+1 + blbl_^_^aiai+i . 


( 4 ) 
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It is important to stress that the Hamiltonian H con¬ 
serves the total number of particles TV, which defines the 
only free parameter of the theory and is conveniently 
expressed as the filling p = N/2L. H features two ad¬ 
ditional, relevant symmetries, namely (i) the subchain- 
parity Px = (— 1 )^* (x = a, b), and (ii) time-reversal 
symmetry P = K represented by complex conjugation 
K and = x\^\ 

Ground states — In order to derive the ground states 
analytically, the observation that Hamiltonian H is the 
sum of local projectors and therefore a locally positive 
operator is crucial. Then we exploit the fact that any 
zero-energy ground state must be annihilated by all local 
terms in (1) and (3) simultaneously. That is, if we find 
a state with zero energy which is annihilated by all local 
terms, we can be sure that it is a ground state. This 
yields a viable method to construct them from scratch 
— provided zero-energy ground states exist. 

For an open ladder, there are exactly two degenerate 
zero-energy ground states for each filling 0 < N < 2L, 
Fig. 2 (a), denoted as \N,a) and characterized by the 
upper chain parity a = Pa G {-l-I,— 1 }, see the supple¬ 
mental material for a rigorous proof. For an appropriate 
fermion gauge, see Fig. 2 (b), each ground state is given 
by the equal-weight superposition of distributing N par¬ 
ticles on the two chains constrained by the fixed subchain 
parity a. To cast this in a formal description, we first in¬ 
troduce the fermion number states \n)^ with x = a, b 
and n G { 0 , 1 }^, i.e., \n)^ = {a\)P.. (aJ)T.. | 0 )„ 

for the upper chain with the number of fermions \n\ = 
J2i=i Then the equal-weight superposition states on 
each chain with a fixed number of particles reduce to 
l^)a: = S|n|=n I”)a:’ State is not normal¬ 

ized. Finally, the equal-weight superposition with fixed 
particle number N and subchain-parity a can be written 
as 

\N, a) = En, (-!)"=« \n)a 1^ - n)b (5) 

where jv a is the normalization factor that counts the 
number of superimposed fermion configurations. 

In contrast, for a closed ladder the situation is more 
subtle: For even total particle number N = 2K, there is 
a unique zero-energy ground state \2K,—1) in the odd- 
odd (a = —I) subchain parity sector, whereas in the 
odd-iV sectors all states are lifted to hnite energy. This 
is summarized in Fig. 2 (c) where the low-energy scaling 
is given as well (see below). 

At this point it seems advisable to compare these 
ground states with those of a single Majorana chain 
(Kitaev’s chain), which in analogy features two zero- 
energy ground states for open boundary conditions [ 8 ]: 
For vanishing chemical potential (perfectly localized edge 
modes), the ground states of the Majorana chain are 
given by the equal-weight superposition of particle num¬ 
ber states with fixed (global) parity. In contrast, here the 


(a) 


p = +1 




(b) Fermion ordering (gauge) 

^-ooooooo<>_; 

r'-oooooooo^ 


Total Parity 

P = +1 _ P = -1 (c) 



Figure 2. Ground states, (a) For every filling N with parity 
P = (—1)^, there are two degenerate zero-energy ground 
states for open boundary conditions, characterized by their 
(upper) subchain parity a = Pa- (b) The chosen fermion 
gauge leads to the simple description of the ground states 
given in the text, (c) Behavior of the spectrum in the low- 
energy sector of symmetry subspaces classified by the total 
parity P and the subchain parity a. Both open (OBC) and 
periodic (PBC) boundary conditions are shown, zero-energy 
states are drawn bold. 


chains act as mutual particle reservoirs and the ground 
state degeneracy arises due to two admissible subchain- 
parity configurations within each fixed particle number 
sector. 

We start exploiting the concise description of the 
ground states, and derive simple expressions for density 
correlations, superfluid order parameter and the Green’s 
function (single particle correlation). To this end, it 
proves useful to define the parity-split binomial coeffi¬ 
cients (PsBC) which count the configurations to dis¬ 
tribute N particles among EEi sites with the ad¬ 
ditional constraint that the parity of subsystem Li (1 < 
j < (?) is fixed by ai = ± 1 , 


( Lu-,L, \ 


N 

E 


Lg 


nEi (6) 


with = [1 + Q!i(—1)"’] /2. Although we are not aware 
of simple analytical expressions (except for special cases, 
see supplement), the PsBCs can easily be evaluated nu¬ 
merically. Due to the simple structure of the ground 
states, all correlation functions and expectation values of 
IA^, a) can be rewritten in terms of finite combinations of 
PsBCs. E.g., the normalization of the two ground states 
reads Nl.n.c = 

We find that the density-density correlation function 
factorizes, {xlxiyjPj) —)■ for i 7 ^ j; x,y G {a,b} in 
the thermodynamic limit L, A^ —> 00 with fixed particle 
density p. The pair correlations read \ {x\x\j^pyjyj+i) \ -G 
p^(I — pY ior i Y j \ tx,y G {a, 5}, and indicate a con¬ 
densate of p-wave pairs with true long-range order. Note 
that the results for both correlators do not depend on the 
subchain parity a of the ground states. This is true up 
to exponential corrections vanishing with L —>■ 00 . For 
particularly symmetric setups (e.g., x Y V N odd) 
these corrections even vanish identically. 
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Figure 3. Ground state properties, (a) Intra-chain single 
particle correlation {a\aj) (Green’s function) as a function of 
the distance \i — j\ for various fillings p and a chain of length 
1/ = 30. The revival for \i — j\ ~ L indicates exponentially 
localized edge states (grey region), (b) Overlap of the ground 
states for time-reversal invariant (TRI) and breaking (TRB) 
perturbations of H in dependence of the position i,j of the 
subchain parity violating single-particle hopping (blue: p = 
0.5 TRI, red: p = 0.5 TRB, grey: p = 0.25 TRB). 

The intra-chain Green’s function (indicating single 
particle off-diagonal long-range order [29]) can be ex¬ 
pressed in terms of PsBCs {j > i -b 1) 


ground state degeneracy for the donble wire setup can 
either be protected by time-reversal symmetry T or sub¬ 
chain parity Px, and is only spoiled if both symmetries 
are broken at the same time. The latter, however, is not 
surprising as the two edge states on the upper and lower 
wire are not spatially separated. We will show below that 
our model can be generalized to wire networks, where the 
different edge states become spatially separated. Then it 
follows immediately that the topological properties are 
protected against any local operator O conserving the 
total number of particles. 

Ground state entanglement — Another well-known sig¬ 
nature of topological states is a stable degeneracy of the 
entanglement spectrum (ES) [31-33]. In our case, the ES 
of the ground states \N,a) with respect to a bipartition 
(S|L \ S) of the ladder [see inset of Fig. 4 (b)] is given by 
the Schmidt decomposition 

1^: a) =51 |n,^)glAf-n,a/3)L\g (10) 

n /3=±1 


—-^L.N.a [A-(-l.-a ^-l,a] (7) 

where n-i supple¬ 

ment for a detailed derivation. In the thermodynamic 
limit one finds exponentially decaying correlations in the 
bulk, see Fig. 3 (a), 

(a;|a;j) = for 1i, j L; x G {a,b} ( 8 ) 

where 7 is some function of the filling with 0 < 7 (p) < 00 
and 7 ( 1 / 2 ) = 00 . The boundary terms read |(a|aL)| —)■ 
p{\ — p) in the thermodynamic limit, indicating the ex¬ 
istence of exponentially localized edge states. Fig. 3 (a). 

The topological protection of the ground state degener¬ 
acy is most conveniently characterized in terms of their 
indistinguishability by any local perturbation [14, 30]. 
Let O be an arbitrary local (hermitian) operator. Then 
the expectation values (aj O |a) and (—a] O |—a) are 
identical up to an exponentially small correction — as fol¬ 
lows from the above analysis of the correlation functions. 
However, for operators violating the subchain-parity P^,, 
also the overlap (—a] O |a) must be taken into account. 
Then the situation is more subtle. We illustrate this for 
the simplest case of a single-particle inter-chain hopping 
(the statements can be generalized to more complex Px- 
violating terms, though). Let Oj = 
with complex hopping phase (j) G [0,27r). Splitting 
this perturbation into time-reversal invariant (TRI) and 
breaking (TRB) contributions, one finds by evaluating 
the corresponding PsBCs 

TRI : (—a] agbg + bgag |a) —> 0 (9a) 

TRB : (—a] iogbs — ib\ag |a) —>■ (9b) 

for the distance 6 from the edges of the ladder, 6 L 
when L —>■ 00 and p is fixed. These site-dependent over¬ 
laps are illustrated in Fig. 3 (b). Thus the topological 


and can be written in terms of PsBCs 


Cn,/3 — 


In 


(L-S,L-S\ ,'S.SN ,(L,L\ 

V a/3 ' N-n^ P ' n' ^ « /N 


( 11 ) 


where max{0,fV — 2L -|- 25'} < n < min{A^, 25} and 
/I = ±1. The /? = ± 1-branches of the spectra for a half- 
split system of length L = 20 are shown in Fig. 4 (a) for 
different fillings N and reveal the two-fold degeneracy of 
the ES due to the subsystem subchain parity (3. 

In addition, the scaling of the entanglement of a sub¬ 
system S with the environmental system as a func¬ 
tion of the subsystem size 5 in terms of the entangle¬ 
ment entropy 5®"* [S] = — Tr [ps In p§] , with reduced den¬ 
sity matrix ps = TrL\s [p], yields insight into the low- 
energy physics of the theory. 5®''* can easily be com¬ 
puted from the entanglement spectrum via 5®"*[S] = 
Y^n /3 Fig- 4 (b) shows the (filling dependent) 

variation of entanglement between a growing subsystem 



Figure 4. Entanglement, (a) Two branches (/3 = ±1: 
red/blue) of the entanglement spectrum for a chain of length 
L = 20 and splitting 5 = 10 with fillings N = 10, 20, 30 (dia¬ 
monds, circles, squares). The half-filling branch is highlighted 
grey. Physically, the index n describes the subsystem filling 
while j3 describes the subsystem subchain parity. This illus¬ 
trates the two-fold degeneracy of the entanglement spectrum, 
(b) The entanglement entropy 5®"* as a function of subsys¬ 
tem size S for various fillings p. It obeys an area law with 
logarithmic corrections. 
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and its environmental system as function of the subsys¬ 
tem size S: It obeys an area law with logarithmic cor¬ 
rections, as expected from a critical (gapless) ID system. 
That is, in contrast to the gapped Majorana chain, here 
we face a low-energy theory of gapless Goldstone modes 
due to particle number conservation. With this in mind, 
we have a closer look at the low-energy excitations. 

Low-energy excitations — The single-chain Hamilto¬ 
nians for an open ladder can be mapped to the 
ferromagnetic, isotropic Heisenberg chain via a Jordan- 
Wigner transformation. The complete spectrum of is 
therefore accessible via the Bethe-Ansatz [34]. Exploit¬ 
ing this mapping, it is possible to construct the analog of 
single magnon states for our theory. These exact low en¬ 
ergy eigenstates for the open double-ch.am take the form 

jfc; TV, a) = 0 IA^, a) (12) 

with momentum k = 0 < m < L, and the operator 

= EU cos [|(2j - 1)] (- 1 )D^“. . (13) 

The eigenenergies are given by a quadratic excitation 
spectrum Ek = 4 sin^ |. This behavior of the Goldstone 
mode is in excellent agreement with the appearance of a 
true condensate and vanishing compressibility; recall that 
for any fixed number of particles there is a zero-energy 
ground state. An equivalent behavior is well-known for 
non-interacting bosons and the ferromagnetic Heisenberg 
model in one-dimension. The interpretation of these fea¬ 
tures is that our model is exactly solvable at a critical 
point. 

Wire networks and non-abelian statistics — A crucial 
aspect of our model is that the derivation of the exact 
zero-energy ground states can be straightforwardly gen¬ 
eralized to much more complicated wire networks con¬ 
sisting of open and closed single chains sectionally con¬ 
nected to ladder segments with arbitrary positive cou¬ 
pling strengths, see Fig. 5 (a) for an example; the gen¬ 
eral formalism is presented in the supplement. It follows 
immediately that the ground state degeneracy scales as 
2 E/ 2-1 E > 2 the number of open subchain ends. 
This scaling is in agreement with the interpretation of 
the localized edge states as interacting equivalent of Ma¬ 
jorana zero modes. In order to provide a rigorous proof 
of the topological properties characterizing the localized 
edge states, we derive the full braiding statistics. Note 
that the gap A closes algebraically, A oc I/T^, Fig. 2 (c). 
This still allows for a generalized notion of braiding and 
thereby probing the edge state statistics [14]. 

In order to braid two localized edge states, we con¬ 
sider the wire network of four open subchains coupled 
by a common “bath” chain depicted in Fig. 5 (a) and 
described by Hq. Only the highlighted chains Li {i = 
1,... ,4) take part in the braiding evolution. Thus the 
greyed out subchain can henceforth be neglected and 







Figure 5. Braiding, (a) Setup of four open chains Li, 
i = 1, 2, 3,4 (black) with controllable weak single-particle cou¬ 
plings between the inner four endpoints. The partner chains 
(grey) are not involved in the braiding and can be disregarded, 
(b) The dynamics takes place in the 8-dimensional Hilbert 
space spanned by the subchain-parity eigenstates with fixed 
total parity a = 01020303 = —1. The colors denote the 
subchain parities Oi of the four black chains, (c) Spectrum 
of the weak coupling Hamiltonian during the braiding proce¬ 
dure depicted below the plot. A black arrow indicates sin¬ 
gle particle hopping. There are four degenerate zero-energy 
ground states. The deviation from zero-energy (perfect adia- 
baticity) due to the finite time evolution is shown in the inset 
(~ 10“®). (d) Time evolution for the initial zero-energy state 
|00;0). Shown are the (moduli of the) overlaps with [00; 0) 
and |11;0) (see inset). 


considered as a “bath”, the effect of which is fully in¬ 
corporated into the exactly known ground states. Note 
that the zero-energy states of the uncoupled subchains 
are given by the total filing N and the subchain parities 
ai,..., 04 , spanning a 2 ^ = 16 dimensional ground state 
space in each particle number sector. As we are only 
considering interactions between the four subchains, the 
total subchain parity a = J[[j Oi is conserved and may 
be hxed at a = — 1 , reducing the number of relevant 
ground states to eight, see Fig. 5 (b). The braiding of 
the edges states is described by and achieved by 

adiabatically turning off the coupling between two edges 
and turning on the coupling between the next two edges; 
the full sequence of couplings for the winding of two 
edge states around each other is shown below Fig. 5 (c), 
where arrows indicate single-particle couplings analogous 
to Ai{\ + Aj). 

The analysis is performed by the full numerical time 
evolution of the Hamiltonian H{t) — Hq 0 e/L^iJint(t) 
with e 1 and 0 < t < 8 to guarantee the (quasi)- 
adiabatic evolution. Starting with the initial zero-energy 
state [00; 0), Fig. 5 (d), characterized by oi = —1 = 
02 and 0304 = — 1 , yields the orthogonal final state 
[11; 0 ) = exp [—i / dt iJint(t)] [ 00 ; 0 ), characterized by 
ai = +1 = 02 and 0304 = —1. Repeating the analy¬ 
sis for alternative braiding operations, we find the non- 
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abelian holonomy acting on the degenerate ground state 
space that qualifies the edge states as Ising anyons [35], 
which corresponds to the braiding statistics of Majorana 
edge modes in non-interacting theories. 

Conclusion — We presented a microscopic model of 
interacting fermions giving rise to a gapless topological 
state with non-abelian edge states. The system is at a 
critical point and certain perturbations to the Hamilto¬ 
nian will drive the system into a phase separated state 
(e.g., increasing the attractive interactions), while we ex¬ 
pect resilience of the topological properties against other 
perturbations (e.g., increasing the hopping). Then the 
ground state should be well described by an approach 
based on bosonization similar to [22-26], and might be 
connected to the state studied with DMRG [27]. 

Note added — During the final steps of preparation, 
we became aware of related results studied by lemini et 
al. [36]. 
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Derivation of zero-energy gronnd states 

It is straightforward to show that the two states |-/V, a), 
a = ±1, in each filling sector N are annihilated by H, 
and therefore belong to the ground state space. Here 
we sketch how to derive their uniqueness by constructing 
the zero-energy eigenstates of H from scratch. To this 
end, consider the binary vectors {n,m) G { 0 , 1 }^^ de¬ 
scribing the fermion configurations on the two-leg ladder. 
The Fock space is spanned by number states \n,m) = 


\n)^ \m)f^ G H, where we choose the fermion ordering 
(gauge) according to Fig. 1 (a); this defines a represen¬ 
tation of the fermion algebra {a^, aj, 6 ^, 6 j}i<i<L on H. 
The action of this representation allows us to write 

(t+A^) \n, m) = \n, m) - \e^(n, m)) ( 1 ) 

and 

+ Bi) |n, m) = |n, m) - |ri(n, m)) (2) 

where we used the single-bit swap-operators 


e“((...,ni,nj+i,...),m) 

for arbitrary (n, m) and the bit-pair swap-operators 

Ti((. . . jTli, ...),(... 

for (n, m) with rii 

It is convenient to introduce the partition = 

{ 0 , 1 }^^ of binary vectors with the decomposition Hat = 
B^ U BJj, where is characterized by vectors {n,m) 
with fixed filling + rrii) = N, and fixed sub¬ 
parity = a. It is now easy to see that (i) 

ef,Ti : B'^ -G Bff are bijections for a: = a,b and 


We proceed with a generic state Ilf) GPL, 


((... 

,ni+i,ni, ...),m) 


(3a) 

(n,(. 

..,TOi+i,mi,...)) 


(3b) 



■■)) 

(4) 


1 < i < L (note that efef = Id = nn), and (ii) for 
any pair (n, m), (n', m') G Bf^ there is a finite sequence 
a of ef,Ti such that [n,m) = a{n',m'). That is, the 
connected components of { 0 , 1 }^^ under the described 
family of bit-operations H = {tt} = are 

exactly the B'^ for N = 0,... ,2L and a = ±1. 


rii+i and rrii = rrii+i, and Ti = Id otherwise. 


1^) = (5) 

and evaluate the energy expectation value 

(d/| H |vk) = E.,. ('I^l Af{t + ) l^k) + E. ('PI + Bi) Ivk) (6) 

which reads 


1^-) = E,rGnE(n,r«),K,r«') m')(n, m) {n',m'\[\n,m) - | 7 r(n,m))] (7a) 

= E^GnE(n,r«) - ^*7r(n,m)5'(n,m)] . (7b) 

Doubling the sum, substituting 7 r(n, m) -G {n, m) and using the bijectivity, 7 r{ 0 ,1}^^ = {0,1}^^, yields 

(^-1 77 1 ^-) = i m)«'(n, m) - ^'* 7 r(n, m)^'(n, m)] ( 8 a) 

+ I - '^i*{n,rn)^n{n,m)] ( 8 b) 

= I E^GnE(r^,m)l^(”."i)-«' 7 r(n,m)|^ . ( 8 c) 
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Splitting the second sum into connected components gives the hnal result, 

J2n,c. Evren "1)1^ . (9) 


Since |'I') is a zero-energy ground state if and only if 
(d>|iJ|'h) = 0, and the are the connected compo¬ 
nents under the action of 11, we conclude that for any 
normalizable zero-energy state it is necessary and suffi¬ 
cient that the amplitudes are constant within each sector 

TDOi 

H I'h) = 0 4'(n, m) = V(n, m) G . (10) 

These are exactly the equal-weight superpositions | A^, a) 
of fixed particle number N and subchain parity a intro¬ 
duced in the main text. We stress that these statements 
are modified for closed boundary conditions within the 
even-even and odd-even/even-odd sectors where some of 
the summands in Eq. (9) read m)-|-tE'7r(n, m)| 

due to the fermionic statistics. From this it follows imme¬ 
diately that no normalizable state with exact zero-energy 
can exist, though the ground state energy may vanish in 
the thermodynamic limit nevertheless (see below). 


Parity-split binomial coefficients 


As explained in the main text, the evaluation of ar¬ 
bitrary expectation values and correlators of the zero- 
energy ground states \N,a) is particularly efficient due 
to their simple structure as equal-weight superposition 
of number states with fixed total particle number N 
and (upper) subchain parity a. The evaluation is then 
reduced to counting fermion configurations with spe¬ 
cific parity constraints Oi,..., ctg-i on g — 1 subsystems 
Li,..., Lg_i (the (?th subsystem is the “rest” or “en¬ 
vironment”). Therefore we introduced the parity-split 
binomial coefficients (PsBC) 


( ) 


N 

E 


ni,...,ng-i 


( 



5-1 



( 11 ) 

with = [1 + ai(—1)"’] /2. They can be evaluated effi¬ 
ciently on a computer, though we are not aware of closed 
forms for arbitrary splittings g. However, the (most im¬ 
portant) PsBC used for the normalizing factor and de¬ 
scribing the splitting of the double chain in upper and 
lower wires of size L, AfL,N,a = allows for a closed 

form, namely 


(“)» = Iff) + f ■ (12) 


Thus the leading order is \ (^). This is exact for odd N, 
as one would expect due to the chain-exchange symmetry 
relating the even-odd and odd-even sectors. For even 


fillings N, however, there are corrections of order (^^ 3 ) 
due to the inequivalence of the even-even and odd-odd 
sectors. 

As a final remark, we give the following alternative 
form of a generic g = 2 split PsBC 


( Li,L2\ 

V a ! N 


(^^+'^=) + fELo(-i)”(t)(/-!„) (13) 


which separates the leading contribution from the a- 
dependent corrections. We make use of this result below. 


Evaluation of correlators 


To exemplify the evaluation of correlators, we show 
how the density-density correlations, the suprafluid pair 
correlations, and the intra-chain Green’s function can be 
evaluated in terms of PsBCs. 

We start with the density-density correlations for 
which we find 


{x\x,y]yj) = 


(L,L 

V —a. } iV—2 ^ ^ 

/ L —2,Z/\ /L,L 
\ a ) N-2 \ a J N 


-1 


-1 

N 


for X ^ y 
for X = y 


This follows from counting all fermion conhgurations in 
\N, a) with occupied sites i and j. Along the same lines, 
counting configurations with occupied sites j and j -f 1 
and vacant sites i and * -I- 1 yields similar expressions for 
the suprafluid pair correlations. 


\i44+iyjyj+i)\ = 


L —2,L —2\ (L,L\ 

a } N—2 \ a / 
L-4,L\ 

a J N-2\ a J N 


-1 

N 


for X ^ y 
for X = y 


In the light of these two correlators, it seems advisable 
to consider the general expression 

/ ij — ll,L — ^ 

V a J ]SI—m\ P /jv 


in more detail. Using Eq. (13) from above and setting 
N = 2Lp, one can derive the expression 


V a l2Lp-m\ P ' 2Lp 


C 2 Lpr :^)&" (14) 

aO{e-^^)+pO{e-f^^) 


which shows that any distinction due to the subparities 
a, /3 vanishes exponentially. Straightforward simplifica¬ 
tions for L —>■ 00 lead to the final result 


^L—li,L — l2\ ^L,L\ 

\ a. / 2Lp—m\ 0 / 


-1 

2Lp 


(15) 


p™(l-p) 




+ aO{e-^^)+ 
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Figure 1. Green’s function, (a) Fermion ordering (gauge) defining the number state basis used to construct the ground 
states, (b) Derivation of the single particle intra-chain Green’s function of the zero-energy ground states |A^, a) by counting 
states of fixed subsector parities ai and ol 2 . The single particle hopping from site f to the (vacant) site i measures the parity 
ai of the subsector in between (yellow), (c) The partition used for the derivation of the inter-chain single-particle overlap 
{N, —a\ albj |Af,a). The single particle hopping between site j and the (vacant) site i measures the double-segment (yellow 
and red) parity 0203 due to the fermion gauge shown in (a). The results in both cases can be expressed in terms of PsBCs 
and are described in the text. 


The correlators from above follow now as special cases, 
{xlx^y]yJ) = and |( 44 -Hi% 1 /i+i)l = up to 

exponential corrections depending on the subchain parity 
a of the ground state \N,a). 

To conclude this section, we focus on the intra-chain 
Green’s function. Here the derivation is more subtle since 
the parity of the single-chain segment separating the two 
support regions of the operators must be taken into ac¬ 
count. For j > i one finds 

{N, a\a\aj |7V, a) = [A+i(16) 

\ ^ J N 

where 


The total parity of the first two segments is fixed to 
aiQ !2 = —a since a single fermion is excluded and used 
for the hopping. Therefore we find exactly A-|_i__a config¬ 
urations which contribute - 1-1 (the parity of the segment 
between i and j is -|- 1 ) and A_i^q, configurations which 
contribute —1. Divided by the total number of configu¬ 
rations in \N,a), this yields the result presented above. 

A numerical evaluation is given in Fig. 3 (a) of the 
main text and reveals the single particle off-diagonal long 
range order close to the boundaries of the ladder due to 
the fixed subchain parity. 

Inter-chain single-particle hopping 


A/3i./32 



l,L-j + i-l,L' 
/3i, ^2 


(17) 


is a g = 3 split PsBC. This follows from simple count¬ 
ing arguments and the convenient partition shown in 
Fig. 1 (b): alaj chooses all states of the equal-weight 
superposition in |A^, a) with vacant site i and occupied 
site j. These states are not annihilated and contribute 
±1 to the result as there is exactly one bra state in {N, a\ 
for each number state in a-Oj |Af, a). The sign of the con¬ 
tribution depends on the parity of the segment between 
site i and j due to our fermion ordering from left to right 
along the chains. Fig. 1 (a). We therefore need to count 
all configurations of the — 1 unaffected fermions (one 
is fixed by the hopping) on 5 = 3 segments: 


1. Segment: length \i — j\ — 1 and parity ai between 
sites i and j. 


2. Segment: the rest of the upper chain of length L — 
\i — j\ — 1 (excluding sites i and j) with parity a 2 . 

3. Segment: the lower chain of length a with fixed 
parity a = a(— 1 )^. 


We extend the previous section by a detailed deriva¬ 
tion of the subchain parity violating single-particle hop¬ 
ping between adjacent chains. As already explained in 
the main text, it is convenient to split the the generic 
perturbation 

Oj=e^'^a]bj+e-^’l>b]aj (18) 

with (j) G [0, 27r) into time-reversal breaking (TRB) and 
time-reversal invariant (TRI) contributions, 

TRI : h+ = a]bj+b]aj (19a) 

TRB : hj = a]bj - b]aj . (19b) 

We are interested in overlaps {—a\h^ |q;) where we omit 
the invariant total particle number in the naming scheme 
of the ground states. This reduces to the computation 
of the inter-chain single-particle correlators (—q;| Ojbj |a) 

and (—aj b^^aj |a). 

Splitting the double chain into four segments, as 
sketched in Fig. 1 (c), allows an evaluation along the 
same lines described above for the intra-chain Green’s 
function: 
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Counting fermion configurations with one of the sites i,j empty and the other occupied, sorting them according to 
their subsegment parities i = 1,2, 3, yields the expression 


(-a| a\bj |a) — ^ ^ {^a) n ^ (^“.+i.+i + A-a,-i,-i - Aa,+i _i - A_q__i_+i) 

with the 5 = 4 split PsBCs 

Aqi, 02,0,3 = ( ’ai,ala3 N-1 ' 

The special case i = j allows us to write 

(—a| |a) oc (Aq;,+i.+i + A_q,,_i,_i — Aq,,+i._i — A_q,,_i,+i) ± (Aq,,_i,_i + A_q,,+i,+i — A_q,,+i,_i — Aq,,_i,+i) 

( 22 ) 

for the Hamiltonian perturbation overlaps. Doing the math yields for the time-reversal invariant (TRI) perturbation 


( 20 ) 




\ 1“/ Z^ni,n2,n3=0 VAT-l-ni-na-ris/V m / V na A 713 M ^ 2 

Incorporating the normalization factors yields 

(-a| hi \a) = [1 - (-1)^] (-1)^ - [1 - (-1)"^] (-1)^ O (e-^W^) 


(23) 


(24) 


for N = 2pL, L —>■ oo where we used expression (12) for the normalization factors. Therefore we find that the overlap 
due to TRI perturbations vanishes identically for even fillings N, whereas in odd filling sectors exponentially vanishing 
corrections occur. 

Finally, an analogous calculation for the time-reversal breaking (TRB) hopping yields 

(-a|).- io) ,xa („_,_;r-™-..)(0(A)(0 (-i)”-”’=» e"- (-i)"(‘b‘)(,;tTii) (25) 


which then gives rise to the lengthy expression 

{-a\h-\a) =2a ^ 


/ L \2 ' / 2 j- 2 \ / 2 L- 2 j \ 

\N/2J Z^n=0 \ \ n ) \N-l-nJ ' 


- 1/2 


(26) 


A numerical evaluation was given in Fig. 3 (b) of the main text and shows for fixed filling p and L —> 00 exponentially 
decaying overlaps in the bulk (the localization length depending on p, vanishing at half-hlling) but finite ones for 
j —>■ I, L. This is derived easily in the odd filling sector where for j = 1,L one Hnds 


I (-a| K.l |a) I = 2 Cn-i) Cn) 
which becomes 1/2 for half hlling, p = 0.5. 

The vanishing overlaps for TRI and the finite, edge- 
localized overlaps for TRB perturbations are analogous 
to the fusion of Majorana fermions in a parallel two- 
wire setup of non-interacting Majorana chains due to 
time-reversal symmetry breaking couplings — which il¬ 
lustrates the breakdown of a Z topological index in sym¬ 
metry class BDI to the Z 2 index of D [1, 2]. 


Low-energy sector 

Here we explain the gap and ground state energy scal¬ 
ings depicted in Fig. 2 (c) of the main text in more 
detail. We start with the algebraic energy gap closing 
AEo{L) (X in all four parity sectors (total parity 


~ 2 p(l — p) for L —)■ 00 (27) 


(—1)^ = ±1, upper subchain parity a = ±1) for open 
boundary conditions. In this case, there is a unique zero- 
energy ground state |A^, a) in each sector and we ask 
about the energy of the first excited state defining the 
energy gap AEo{L). As explained in the main text, one 
can derive the single particle excitations in each sector 
{N, a) rigorously. The lowest lying single-particle excita¬ 
tion has momentum = J and energy 

£;i(L)=4sin2^^^ (28) 

which yields an upper bound for the energy gap, 

AEo{L)<^. (29) 
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Figure 2. DMRG Results, (a) Log-log plot of the energy 
gap AEo scaling for L —>■ oo for open and closed ladders in 
the even particle number sector within the odd-odd subsector 
for p = 0.5. Comparison with the analytic result for single- 
magnon excitations with momentum k = ix jL (OBC) and k = 
27r/L (PBC) shows perfect agreement. The insets show the 
same data in linear scale. Note that for OBC the results can 
be transferred to the other subchain parity sectors, (b) Log- 
log plot of the ground state energy Eq scaling for L ^ oo for a 
closed ladder in the odd particle number sector for both fixed 
filling p = 0.5,0.25 and particle number N = 5. For fixed 
filling, the ground state energy vanishes with 1/L, for fixed 
particle number with l/l?. Note that exact diagonalization 
(ED) up to 1/ = 15 suggests a finite ground state energy 
(and thus a single particle gap) which is then invalidated by 
DMRG up to L = 108. The simulations were performed with 
the ALPS libraries [3, 4]. 


Note that we did not disprove the existence of potentially 
lower lying bound states. 

This is, in general, a highly nontrivial task. For in¬ 
stance, it has been shown rigorously that the gap above 
the zero-energy ground state for the isotropic Heisen¬ 
berg chain is given by 1 — cos(7r/L) = 2sin^7r/2L [5], 
which indeed agrees with the single particle excitations 
for A: = ^ given above (the additional factor of 2 fol¬ 
lows from the definition of our Hamiltonian in terms of 
spinless fermions). From a rigorous point of view, this 
statement can only be transferred to -I- which de¬ 
scribes two non-interacting isotropic Heisenberg chains. 


However, our results suggest that the additional interac¬ 
tion H°‘^ cannot decrease the first excitation energy given 
by the k = tt/L magnons (which would certainly be true 
for weak perturbations e ^ 1). This statement 

is underpinned by DMRG simulations for half filling up 
to L = 108 with very high accuracy, see Fig. 2 (a). We 
conclude that for OBC the energy gap closes as 

AEo{L) ~ ^ (30) 

in all sectors {N,a). 

We proceed with the case of periodic boundary con¬ 
ditions. As explained in the main text, there is a dif¬ 
ference between the odd-odd sector ((—1)^ = -|-1 and 
a = — 1) with an exact zero-energy state and the remain¬ 
ing even-even and even-odd/odd-even sectors with finite 
energy ground states. Consider the odd-odd sector first. 
There the fermionic statistics has no effect and the en¬ 
ergy expectation value has the form of Eq. (9) for OBC. 
Performing a Jordan-Wigner transformation yields the 
isotropic Heisenberg chain with (untwisted) PBC. Due to 
the survival of the SU(2)-symmetry, all statements about 
the ground state and the single-particle excitations carry 
over from the previous case of OBC; in particular the gap 
closing, now with AEq{L) ~ 47r^/L^ since the lowest- 
energy magnon has momentum k = 211 jL (cf. fc = tt/L 
for OBC). These statements can again be verified to high 
precision with DMRG simulations, see Fig. 2 (a). 

We conclude with the sectors without zero-energy 
ground state, i.e., even-even and even-odd/odd-even. 
Here we aim at the scaling of the ground state energy 
Eo{L) oc L~^ compared to the zero-energy ground states 
in the odd-odd sector. In contrast to all previous cases, 
here we have neither access to exact ground nor excited 
states which complicates the analysis considerably. 

We start with the even-even sector. Here one can show 

§ > Eo{L) > g (31) 

rigorously. The upper bound follows with the ansatz 
wave function {0 < K < L) 


I Go) = exp 


. TT 



L 

^s(ajas 

S = 1 


bibs) 


\ 2 K,+ 1 ) 


(32) 


which imprints a single-particle phase field onto the 
equal-weight superposition and thereby satisfies the 
twisted boundary conditions with zero-energy in the ther¬ 
modynamic limit. It can be shown by straightforward 
calculations that 

(Go|i/|Go) = 2.4p(l-p).L.sin2^ 

~ 2p(l - p)7r2 • L-i (33) 

with the filling p — N/2L = KjL. The energy is due to 
intra-chain interactions alone, the inter-chain contribu¬ 
tions due to vanish. 
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Finding the lower bound (thereby establishing the al¬ 
gebraic decay) is much more subtle as it makes state¬ 
ments about all possible states — in contrast to the up¬ 
per bound which was derived by a single ansatz wave 
function. The basic line of thought is to consider the 
energy expectation value in Eq. (9) for a generic state 

(^-l H 1^-) oc Efo-)6G 1^* - (34) 

where now = ±1 due to the fermionic statistics 
and i,j are indices encoding the fermion configurations 
{n,m). G can be thought of as a graph describing the 
possible local transformations 'I'i i—>■ due to the Hamil¬ 
tonian H. One can now construct a set C of length-L cy¬ 
cles C in G such that (i) every configuration i is visited 
by at least one of the cycles and (ii) every cycle contains 
one edge (ij) with = —1, preventing zero-energy on 
each cycle. Physically, these cycles can be interpreted as 
virtual tunneling of a single fermion around one leg of 
the ladder: 


an extensively growing energy expectation due to the 
pair-phase locking of the interaction H°‘^. 

In this picture of single particle-phases (naively sup¬ 
ported by bosonization) one might therefore conclude 
that the ground state energy is finite, limL_>oo Eq{L) > 0, 
which then established a single particle gap, i.e. a finite 
difference of ground state energies between even and odd 
total filling. This conclusion is even supported by ex¬ 
act diagonalization for p — 0.5 up to L = 15 where the 
ground state energy seems to approach a finite gap. How¬ 
ever, DMRG up to L = 108 proves this conclusion wrong 
and suggests clearly a vanishing of Eq{L) with L~^ for 
any fixed filling, see Fig. 2 (b). The non-existence of a 
single-particle gap is presumably owed to the “double” 
criticality of the very peculiar point H under considera¬ 
tion, in that it is not only a gapless phase but also on the 
brink of phase separation. 

Wire Networks 


mH\^) 

^ EcgC E(ij)GC (35) 

An application of Holder’s inequality followed by the tri¬ 
angle inequality yields 

(36) 

where the single aij = — 1 is crucial. Due to the construc¬ 
tion of C based on cluster-conhgurations of fermions, the 
double sum can be recast as sum over all configurations 
i weighted by their number of fermion clusters A^c(*), 

^c(*)|'f'.|' = ^>^. (37) 


In the last step we bluntly bounded the ground state ex¬ 
pectation of clusters by one from below which establishes 
the lower bound. 

The stronger statement Eo{L) oc L~^ is again sup¬ 
ported by numerical simulations and analytic results 
known for the Heisenberg model with twisted boundary 
conditions [6-8] . We stress that numerics suggests in the 
ground state (W) L providing the numerically ob¬ 
served lower bound L~^. Note that a bounded cluster 
number, (W) ^ const, would indicate phase separation. 

Finally, consider the closed ladder with odd total fill¬ 
ing, i.e., the even-odd/odd-even sectors. We point out 
that the previous sketch providing the lower bound L~^ 
only relies on the existence of a single chain with even 
Filing. It is therefore still applicable in the present 
case and restricts any decay of the ground state en¬ 
ergy to be algebraic. However, the former ansatz wave 
function (32) yields a finite energy expectation as only 
one of the two chains requires a single particle phase 
due to the even filling. The naive ansatz [Gg) = 

exp |2Ar -I-1, -l-I) (here for a = -1-1) with 

a single-particle phase only in the upper chain fails with 


For the braiding procedure described in the main text 
(Fig. 5) a much more complicated setup than the “sim¬ 
ple” two-leg ladder was used. Here we comment on this in 
more detail: We show that a much wider class of Hamil¬ 
tonians, defined on “wire networks”, can be described in 
exactly the same way as the simple two-leg ladder; in¬ 
cluding the degenerate zero-energy ground states given 
as equal-weight superpositions and the simple evaluation 
of correlators/observables in terms of PsBCs. 

Examples of possible network topologies accessible by 
these techniques are shown in Fig. 3. In (a) two wires 
coupled by a common closed bath chain give rise to a 
two-dimensional ground state degeneracy, characterized 
by the parities of the two subchains. In (b) a more com¬ 
plex patch from a stacked wire setup is shown. The more 
intricate setup in (c) will be used in the following to il¬ 
lustrate the construction of a generic wire Hamiltonian 
derived from the two-leg version of the main text. 

In fact, for any network of g open and g' closed (single) 
chains, sectionally connected as double wires, one can 
write down a generalized version of the wire Hamiltonian 
E[ such that in any filling sector N^ there are exactly 
2®“^ {g > 1) zero-energy states \N-, ai,..., Og-i) given as 
the equal-weight superpositions of fermion configurations 
with the ith subchain parity fixed at = ±1. In the 
special case oi g = 1, the ground state |A^) for each filling 
N is uniquely determined. The parity of closed single 
chains is fixed at a = —1 to make the ground state energy 
vanish identically, so only open chains contribute to the 
ground state degeneracy. Note that g = G/2 where E 
denotes the number of open subchain ends, as used in 
the main text. 

For instance, consider the g = 2 wire network in 
Fig. 3 (a) where the upper “bath” chain is closed and 
therefore has fixed parity ag = ~1- Then the two open 
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chain parities a\ and are constrained by the total par¬ 
ity via a\ai = —(—1)^. Hence the ground state space is 
two-fold degenerate for fixed filling N^ as illustrated by 
the two possible subparity sectors. Although the degen¬ 
eracy is the same as for the simple ladder setup used in 
the main text, here the edge states are spatially separated 
and therefore resilient even against TRB perturbations. 

To write down the local Hamiltonian featuring the 
properties sketched above, we start with an arbitrary 
wire network £ of fermionic sites (for the sake of clar¬ 
ity, think of the example in Fig. 3 (c)). Formally, 
£ = (V, C,7^) is defined on a (so far unordered) set V 
of vertices s (fermionic sites) as a collection C of unori¬ 
ented chains C (open or closed) and a collection V of 
pairing sections P describing the sectional attachment of 
two chains Ci and C 2 to form a two-leg ladder in some 
segment. 

Note that a single chain C can be oriented in two dif¬ 
ferent ways while preserving the locality of the sites; this 
defines a direction of each single chain (which is a com¬ 
pletely unphysical, that is, gauge choice). In Fig. 3 (c) a 
possible choice is given by the black arrows. Note that 
the only relevant manifestation of this (gauge) directions 
is defined on pairing segments P € P as a binary function 
X ■ 'P ^ via x{P) = +1 (~1) if the orienta¬ 

tions of the two bonded chains are parallel (antiparallel). 
In the following, (abstract) wire networks £ are always 
thought of being augmented by an orientation function 
X- 

After these formal preliminaries, we can write down 
the Hamiltonian that governs the physics on the wire 
network £ as follows: 

Hm = Ecec H'[C] + Epcp H"[P] (38) 

where H' describes the intra-chain physics, 

H'[C]=Y.e^cAe{t+A,) (39) 

with the single particle hopping Ag = Ogaf -f apa\ on 
edge e = (s,p), and H" is responsible for the inter-chain 
interaction. 


H"[P]=Y.f^pBf{x{P) + Bf) (40) 

with the known pair hopping between the two chains 
Bf = alajjQraq + ala^^asUp. Here / = {s,p;q,r) defines 
a “face” on the ladder segment with corners s,p on the 
“upper” and q,r on the “lower” chain. Fig. 3 (c). Note 
that the choice of chain orientations affects only the defi¬ 
nition of H” which is rooted in the invariance of Ag with 
respect to orientation inversion on a single chain (s p), 
in contrast to the variance Bf o —Bf induced by orien¬ 
tation inversion on one of the two chains (e.g., s -O- p). 

We point out that in the framework of the fermion al¬ 
gebra {osjsgv, the gauge transformation Og —>■ iog for 
s G C inverts all signs of xiP) for pairing sections P in¬ 
cluding chain C. In this sense, the orientation arrows of 
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Figure 3. Wire Networks, (a) Simple wire network with 
two open chains connected by a common closed “bath” chain 
fixed at odd parity. The ground states are characterized by 
the subchain parities, (b) Patch of a more complex network 
of stacked chains. The degeneracy grows with the number of 
open chain ends E. (c) Example of a more intricate ladder 
setup with a closed “bath” chain fixed at ao = — 1 subchain 
parity and the two open subchain parities ai, 0 : 2 . There are 
two exact zero-energy ground states in each particle number 
sector characterized by —0102 = (—1)^. This scheme is used 
to introduce a generic version of the two-leg ladder Hamilto¬ 
nian introduced in the main text. 


single chains can be inverted by the U(l)-rotation of the 
fermions on the chain in question. Then it is obvious that 
the choice of orientations (and thereby x) h&s no physi¬ 
cal significance as the resulting Hamiltonians are related 
by gauge transformations. It is important to stress that 
while the choice of orientations is unrestricted, the rela¬ 
tive orientations % are implicitly given by this choice and 
the topology of the wire network. That is to say, there are 
different gauge equivalence classes of Hamiltonians for a 
given wire network £, and only the representatives i7[£] 
(of a specihc class) defined above give rise to degenerate 
zero-energy states. 

For instance, there are network topologies where it is 
not possible to gauge all relative orientations to x{P) = 
-|-I (which would be the naive generalization of our 
paradigmatic two-leg ladder to networks). One such ex¬ 
ample is given in Fig. 3 (c) where the colors of inter-chain 
pair interactions encode the sign of x{P) with respect to 
the orientations given by the black arrows. Note that 
there is no choice of orientations so that all ladder seg- 





















ments feature parallel orientations. 

It remains to be shown that the ground state space 
of H[C] features degenerate zero-energy states given by 
equal-weighted superpositions and characterized by the 
parities of the open subchains in C. To this end, a rep¬ 
resentation of the fermion algebra on the Hilbert space 
is required. This is achieved by defining a number state 
basis of the Fock space, based on a total ordering of the 
fermion sites. A convenient choice for the construction of 
the ground states is to order fermions along single chains 
and parallel to the previously chosen orientations. The 
sequence of the single chains relative to each other is 
not important as long as only parity conserving opera¬ 
tions between the chains are concerned (which is true for 
H[C]). Note that the sequence of chains is important 
for the evaluation of subchain parity violating observ¬ 
ables/correlators such as inter-chain Green’s functions 
and single particle hopping. A possible ordering is high¬ 
lighted in Fig. 3 (c). 

The construction of the ground states follows from 
exactly the same reasoning as described above for the 
simple two-leg ladder, generalized to multiple connected 
chains. The single particle hopping within the chains can 
still be described on the previously defined Fock basis as 

Ae{t + A^)\n) = \n)-\£f,{n)) (41) 

with, e = (i,i + 1), 

Eei- . .,ni,ni+i,. ..) = (.. .,ni+i,ni,...) (42) 

where we no longer split the fermion configuration into 
two subchain configurations n and m but consider the 
total configurations n G {0, with the total (single) 
chain length Ltot = of 5 + 5^ chains. For the 

two-leg ladder one has g = 2, g' = 0 and Ltot = L -|- L = 
2L. 

For the inter-chain pair interaction one finds 

BfixiP) + Bf) \n) = \n) - |T/(n)) (43) 

with 

Tj (. . . ,77-2, , . - ■ , Tlj , 5 ■ ■ ■ ) 

= (44) 

if rii = rii+i and Uj = rij+i, and ry = Id otherwise. 
Now it becomes clear that the sign x{P) in H” compen¬ 
sates for the signs due to fermion pair hopping, which 
are different for hopping between chains of parallel and 
antiparallel order (the latter includes an exchange of the 
pair fermions). The effect of H[C] (depending on the 
chain orientations through x) on the used Fock basis \n) 
(also depending on the chain orientations) is therefore 
invariant under the choice of chain orientations and has 
the same form as for the two-leg ladder. 

In conclusion, the energy expectation still has the 
generic form 

(4-1 H[C] I4-) cx I'f'H - (45) 



Weak couplings (c) 


I : 
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Figure 4. Braiding scheme, (a) Braiding setup with four 
chains Li, i = 1,2, 3,4. The “bath” chain is omitted as it is 
not relevant for the low-energy dynamics. The chosen fermion 
ordering defining the number state basis is highlighted red. Oi 
denote the subchain parities fixed to 01020304 = —1. (b) 8 
basis states of the ground state space in which the braiding dy¬ 
namic takes place. The linear combinations |00; 0) and |11; 0) 
are exchanged by the braiding operation, (c) Weak single¬ 
particle couplings used for the braiding. Blue (red) edges 
denote time-reversal invariant (breaking) single-particle hop¬ 
ping. 


where {ri,n') encodes the possibility to transform the 
fermion configurations n n' into each other by an 
elementary process of H[C] (i.e., a single particle hopping 
within a chain or a pair hopping between chains in a 
bonding region). Note that Eq. (45) is only valid for 
number states with odd subchain parity ac = — 1 for 
closed chains C. 

The zero-energy ground states of the positive semidefi- 
nite Hamiltonian H[C] are therefore given by states with 
constant weights on all connected components of G, that 
is, on fermion configurations with fixed total particle 
number N and fixed subchain parities ai, i = . ,g—l. 

This gives rise to the claimed 2®“^ fold degeneracy in 
each filling sector N, spanned by the ground states 
\N-,ai,... ,ag-i). It is now easy to see that arbitrary 
correlators/observables still can be expressed in terms of 
PsBCs; the splitting degree of which may be much larger 
though. 

A detailed analysis of such wire networks (e.g., effects 
due to perturbations of inter-chain single-particle hop¬ 
ping) might be a promising field for future work. 


Braiding 

The setup used in the main text, see Fig. 5 (a), is an ex¬ 
ample of a wire network, the ground states of which are 
given by the equal-weight superpositions of all fermion 
configurations with total number N and subchain par¬ 
ities ai, i = 1,2, 3,4. The parity of the “bath” chain 
connecting the four open chains can be fixed so that 
Q!ia 2 Q; 3 Ck 4 = —1 or ao = and is of no interest in 
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the following. Due to the previous constraint on the sub¬ 
chain parities, the ground state space can be restricted 

to the eight states ^^^ 2 ) fixed particle number N 

and bath parity oq, see Fig. 4 (b). 

For the braiding we consider the time-dependent 
Hamiltonian 




(46) 


where the time-independent part i^o defines the wire 
network without couplings between the endpoints of 
subchains Li (including the bath chain Lq absent in 
Fig. 4 (a)) and the weak single-particle coupling i?int(i) 
describes the time-dependent braiding procedure within 
the above described ground state sector of Hq (Fig. 4 (c)). 


The Hamiltonian i/int(0 couples the four chains Li, i = 1,2, 3,4, via single-particle hopping as follows: 


Hintit) = 



for 

0 < t < 1 

and 

T = t — 

0 


for 

1 < t < 2 

and 

T = t — 

1 

)Hl^,+rHlll, 

for 

2 < f < 3 

and 

r = t — 

2 

)HI^,+tHI^, 

for 

3 < t < 4 

and 

T = t — 

3 


(47) 


and we set iJint(i) = mod 4) to allow for a double exchange (that is, full braiding) of the two edge states 

located on the inner endpoints si and S 2 of chains Li and L 2 for 0 < f < 8. The chain coupling Hamiltonians are 
defined analogous to the intra-chain coupling Ai{t -\- Ai), namely 

Hfnt = + 0^4 ala + ”^* 3(1 - + n*, (1 - Us^) (48a) 

Hint = + ns^{t - UsJ + - n^J (48b) 

Hint = a^i 0^2 -f a „2 -f n*, (1 - ) -f (1 - n* J (48c) 

HlH = 0*10^3 -f 0^3 0^3 -f nsi(ll - n^a) + na3(ll - (48d) 

where Si denotes the terminating fermion sites on chain Li, see Fig. 4 (a). Note that the single-particle hopping 
in iJb). breaks time-reversal symmetry — which is necessary to create a finite overlap between ground states in the 
thermodynamic limit for couplings between chains with reversed fermion ordering. 

In the ground state basis, the coupling Hamiltonian matrix elements read as follows (here exemplarily for i7^{) 


134/32 


H; 


34 


dittaV _ 
0 . 40.21 


^OL,j4-/^ct ,/3i^Q2,/32—/3a^a4, —/34 ^ 3^14 ( ^ 1 : 10 : 2 ) 


Fi, T2, -^3 — 1, L4 — 1, Lq 
ai, a 2 , o^, —04 

4 


N-l 


Li, L 2 , L 3 — 1 , L 4 — 1 , Lq 

Ol, 02, —03, 04 


N-l 


with the shorthand notation Sa^/S = rii=i ^ai,Pi and the normalizing factor 

Li, L 2 , L 3 , L 4 , Lq 


Ol,02, 03, 04 


N 


Li,L 2 , L 3 , L 4 , Lq 
Ol, 02, 0:3, 0:4 


-1 


N 


P(l-P) 


N-l 


for Li 


Na = 

Using that 

1 . 

Li,L 2 , L^ — 1 , L 4 — 1, To ' 

f3i,l32, 333 , /34 

00 up to exponentially small o, /1-dependent terms, leads to the expressions 

^ " 3a,(4 ‘ 1 ,/3i,/32^aa, —/3a^Q4, —/34 ' ( 0 : 30 : 2 ) 

3a,f4 ‘ 1 —/3i:la2,/^2^c^a,/3a^CK4, —/34 ’ ( ^ 0402 ) 

3a,f4 ‘ t T —/3i:1 q;2 /32,/3a^0^4,/34 ’ (~fOi 02 ) 
3a,(4 ' 1 “t“ —/3i :1 q;2,/32 /3a^0^4,/34 ’ ( ^) * 


/h/h 

/I 4 /I 2 

^1^3 

/34/32 

/h/h 

/34/32 

^1^3 

^4/32 


H, 


Oi 03 \ 

O4O2/ 

Oi03\ 

O4O2/ 

Oi03\ 

O4O2/ 

Oi 03 \ 

O4O2/ 


(49) 

(50) 

(51) 

(52) 

(53a) 

(53b) 

(53c) 

(53d) 


These results are valid for —/> oo for all chains i and fixed filling p. We omit the global factor of 2p(l — p) which is 

justified as it is independent of the parities Oi up to exponentially small corrections. 


We integrate the propagator 


numerically and thereby end up with a non-trivial holon- 


U {t) = exp 


-i / dr Hint (t) 
Jo 


(54) 
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omy Us^s 2 — U{4:) on the four-fold degenerate ground 
state space describing the action of exchanging (i/sisa) 
and braiding the two localized edge states on si 

and 52 - 

The ground state space at t = 0 mod 4 is charac¬ 
terized by a coupling of chains L 3 and L 4 and there¬ 
fore spanned by the four states |ai, 02 ) |a 3 a 4 )_|_ with 
aia 2 a 3 a 4 = —1 defined as 


|ai, 02 ) |Q;3a4)^ = 


aia3\ 

0.40.21 


0 \ —Q!3\ 
—04 02 / 


(55) 


up to normalizations. The action of the braiding is 
then diagonal in the 0-404 blocks and reads 

Wi, 02 ) 10304 )+ = 1 - 01 ,- 02 ) 10304 )+ (56) 


Pairing the Majorana modes of two different Majorana 
chains into a fermion mode, one finds that the symmetric 
(antisymmetric) linear combination of the two parity sub¬ 
sectors admissible by the fixed total parity corresponds 
to an empty (occupied) mode. Translated into our case 
this means that the occupation m of the “fermion mode” 
created from “Majorana modes” at si and S 2 is given via 

|m)|m',n34) (58) 

= [|ai, 02) + (- 1 )”" \-oi, -02)] |a 3 a 4 )+ 

which, of course, is just a basis transformation. Here 
0102 = —(—1)™ describes the occupation m' of the 
“fermion mode” defined by pi and P 2 - In this basis, the 
action of the braiding is particularly simple, that is, di¬ 
agonal 

Us,s 2 1”^) N':«34) = (-1)™ \m) \m',n34) (59) 


up to additional phases. 

In order to understand this result more thoroughly, 
it is useful to cast these expressions into the language of 
the original (quadratic) Majorana chain, keeping in mind 
that due to the interactions there are no “fermion modes” 
or “Majorana modes” describing the ground states — 
this is just semantic sugar. 

Remember that the delocalized, fermionic edge mode 
of the original Majorana chain is empty (occupied) if 
the total chain parity is odd (even). In this linguis¬ 
tic framework, we may call the edge state of a sub¬ 
chain in our interacting model empty (occupied) if its 
subchain parity is odd (even). That is, the occupa¬ 
tion ni of the “fermion mode” defined by the pair of 
“Majorana modes” at pi and si, see Fig. 4 (a), corre¬ 
sponds to the subchain parity oi via oi = —(—1)"^; and 
the same for L 2 - This defines the new naming scheme 
I’^ 1 j« 2 ) 1 ^ 34 )+ = |ai, 02 ) 10304 )+ with ni ^2 =0,1 and 
0304 = — (—l)"'^^. Then the action of braiding the “Ma- 
joranas” si and S 2 reads 

C^L.|0^0)|0)+ = |1^1)|0)+ (57) 

which would be interpreted in the quadratic mean-field 
theory of the Majorana chain as creation of a fermion 
pair out of the superconducting condensate by braiding 
two Majorana modes. 


up to a phase independent of m. This is exactly the 
relative phase one would expect from rotating a fermion 
mode by 2 tt with occupation m. 


* nicolai@itp3.uni-stuttgart.de 
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